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We investigate quantum transport through strongly disordered barriers, made of a material with 
exceptionally high resistivity that behaves as an Anderson insulator or a "bad metal" in the bulk, 
by analyzing the distribution of Landauer transmission eigenvalues for a junction where such bar- 
rier is attached to two clean metallic leads. We find that scaling of the transmission eigenvalue 
distribution with the junction thickness (starting from the single interface limit) always predicts 
a non-zero probability to find high transmission channels even in relatively thick barriers. Using 
this distribution, we compute the zero frequency shot noise power (as well as its sample-to-sample 
fluctuations) and demonstrate how it provides a single number characterization of non-trivial trans- 
mission properties of different types of disordered barriers. The appearance of open conducting 
channels, whose transmission eigenvalue is close to one, and corresponding violent mesoscopic fluc- 
tuations of transport quantities explain at least some of the peculiar zero-bias anomalies in the 
Anderson- insulator/superconductor junctions observed in recent experiments [Phys. Rev. B 61, 
13037 (2000)]. Our findings are also relevant for the understanding of the role of defects that can 
undermine quality of thin tunnel barriers made of conventional band-insulators. 
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I. INTRODUCTION 

The advent of mesoscopic quantum physics^^ in the 
early 1980s has profoundly influenced our understanding 
of transport in solids. Advances in micro- and nano- 
fabrication technology have brought about small enough 
structures (< 1/im) in which, at low enough tempera- 
tures (T <C IK), propagation of an electron is described 
by a single wave function (instead of density matrices 
in macroscopic solids) since inelastic dephasing processes 
can be suppressed below the temperature-dependent de- 
phasing length L^. Thus, their transport properties have 
to be analyzed in terms of quantities that take into ac- 
count non-local features of quantum dynamics (such as 
the quantum corrections^ to the conductivity which are 
non-local on the scale of L^), finite-size of the sample, 
boundaries, and measurement set-up of macroscopic ex- 
ternal circuit, rather than using traditional local and self- 
averaging quantities (such as the conductivity) applicable 
to bulk materials at high enough temperatures. 

Particularly influential ideas have emanated from the 
Landauer-Biittiker approacfaJi&l to quantum transport 
which treats conduction within the phase-coherent sam- 
ple as a complicated (multichannel) quantum-mechanical 
scattering problem. This viewpoint introduces a set of 
transmission coefficients as the fundamental property of 
a mesoscopic conductor. The transmission coefficients 
T n are formally defined as the eigenvalues of tt', which 
is the product of a transmission matrix t and its Her- 
mitian conjugate t^. In the two-probe geometry, where 
mesoscopic sample is attached to two semi-infinite ideal 
metallic leads, the t-matrix connects the transmission 
amplitudes of the flux-normalized states in the left lead 
to the outgoing states in the right lead. Thus, the basis of 
eigenchannels, which diagonalizes the matrix ttt, offers 
a simple intuitive picture where conductor can be viewed 



as a parallel circuit of independent transmission channels 
characterized by channel-dependent transmission proba- 
bility T n . Within this framework, pure tunnel barrier 
is a rather simple case where all transmission eigenval- 
ues T„ <C 1 are the same and much smaller than one (the 
opposite limit, T n = 1, is a property of the ballistic trans- 
port occurring through fully open conducting channels). 

Since many electronic devices employ quantum- 
mechanical tunneling through an insulating barrier, their 
design and optimization requires to understand whether 
transport occurs via pure tunneling or if it is affected 
also by the defects in the barrier^S In particular, high- 
critical current density for Josephson tunnel junctions^ or 
impedance level for magnetic tunnel junctions^ require 
ultrathin and highly transparent barriers that can eas- 
ily be pushed out of the genuine tunneling regime. 5 The 
diagnostics of non-trivial barrier properties requires to 
investigate quantities beyond just the conductance since 
its exponential decrease with the barrier thickness, as a 
naive criterion of pure tunneling, can be generated by 
vastly different underlying microscopic mechanisms. For 
example, recent experiments 5 ^ have pointed out how ho- 
mogeneous ultrathin (e.g., thickness ~ lnm) aluminum 
oxide barriers can accommodate high transmission chan- 
nels T n ~ 1 (which are detrimental for various device 
operation 6 ). This is due to extended states induced by 
disorder— or intrinsic transport mechanism in disordered 
mesoscopic systems^ rather than due to rare defects such 
as pinholes with more than unit-cell dimension. 

When static disorder becomes strong enough, solids 
undergo localization-delocalization (LD) transition lead- 
ing to an Anderson insulator^ Such phase is substan- 
tially different from the conventional Bloch- Wilson band 
insulator since density of states at the Fermi energy re- 
mains finite in Anderson insulators. On the other hand, 
the wave function associated with the localized states is 



confined within the region of a characteristic size speci- 
fied by the localization length £. 

Here we explore quantum transport through a strongly 
disordered barrier, separating the two clean metallic elec- 
trodes, by computing statistical properties of the trans- 
mission eigenvalues for an ensemble of three-dimensional 
(3D) samples with different impurity configuration. We 
focus on the appearance of completely open transmission 
channels 11 * 1 ^ T„ ~ 1, as the barrier thickness increases 
from the single interface limit to the junction thickness 
where tunneling through the Anderson insulator takes 
place, and their effect on experimentally accessible trans- 
port properties. That is, the full statistics of T n allows 
us to obtain frequently measured quantities that contain 
the signatures of such non-trivial transparency proper- 
ties: (a) the zero-frequency power spectrum of the shot 
noise^ and (b) the conductance^ Gns of a hybrid junc- 
tion composed of a thin Anderson insulator attached to 
a superconductor, whose unusual properties have been 
unearthed in recent mesoscopic transport experiments. 11 
Our findings on quantum transmissivity of single inter- 
face and thin barriers of a strongly disordered materials 
are relevant also for the analogous classical coherent scat- 
tering problems, such as the light propagation through 
thin, but strongly diffusive, medium^ 

The paper is organized as follows. In Sec. [H] we in- 
troduce the Hamiltonian model of the disordered barrier 
and corresponding real-space Green function technique 
that allows us to obtain an exact transmission matrix 
of a specific sample. In Sec. Mil wc study the scaling of 
the distribution of T n as a function of the barrier thick- 
ness, where disorder strength serves as a parameter whose 
tuning induces Anderson insulator, as well as a "bad 
metal" regime upon approaching the LD transition from 
the metallic side. Measurable transport quantities — shot 
noise and proximity conductance Gns — determined by 
these distributions are discussed in Sec. IIVI In particu- 
lar, we find the shot noise to be a sensitive single param- 
eter characterization of the transparency of multichannel 
barriers, as well as of different types of diffusion through 
dirty metallic barriers. We conclude in Sec. 

II. TRANSMISSION PROPERTIES OF THE 
ANDERSON MODEL FOR DISORDERED 
BARRIER 

We model non-interacting electrons in the disordered 
barrier by a standard Anderson model^S 

H = J2e m \m)(m\+t £ |m)(n|, (1) 

m <m,n) 

which is a tight-binding Hamiltonian (TBH) defined on 
a simple cubic lattice L x L y x L z . The nearest neigh- 
bor hopping matrix element, between s-orbitals (r|m) = 
ijj(r — m) on adjacent atoms located at sites m of the lat- 
tice, is denoted by t and sets the unit of energy. Here L 
is the thickness of the junction in the units of the lattice 
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FIG. 1: The boundaries of different transport regimes, deter- 
mined by the strength of the disorder W , in a bulk 3D con- 
ductor described by the half-filled (Ef — 0) Anderson model. 
At W » 6, the Boltzmann equation breaks down (putative 
semiclassical mean free path becomes smaller than the lattice 
spacing I < a), while at W as 16.5 wave functions become lo- 
calized. Within the intermediate "bad metal" regime, particle 
motion is "intrinsic" diffusion that requires non-perturbative 
quantum descriptions- 



spacing a (i.e., L is equal to the number of disordered 
interfaces of the cross section L y x L z which are stacked 
along the z-axis, chosen as the direction of transport, 
and coupled via hopping t to form the barrier). We set 
Ly — Lz — 20, which yields the quantum point contact 
conductance Gq PC {E f = 0) = 259G Q (Gq = 2e 2 /h 
is the conductance quantum) of the corresponding clean 
system attached to two leads of the same cross section 
(i.e., for this set-up, there are at most 259 fully open Lan- 
dauer conducting channels T n = 1 at half-filling, out of 
400 supported by its cross section^). The disorder is in- 
troduced by setting a random on-site potential such that 
e m is uniformly distributed in the interval [-W/2, W/2]. 
The whole band of the Anderson model becomes local- 
ized, i.e., the LD transition takes place at the Fermi en- 
ergy Ep = of the half-filled band, when critical disorder 
strength W c w 16.5 is reached. 

It is important for subsequent discussion to recall that 
there arc three fundamentally different transport regimes 
in bulk 3D disordered conductors^ (i.e., in the cubes 
L x L x L with a given concentration of impurities): 
(a) the semiclassical regime, where the Bloch-Boltzmann 
theory and perturbative quantum corrections (obtained 
from the Kubo formula) describe resistivity of diffusive 
(I -C L) systems; (b) the "bad metal" regime character- 
ized by exceptionally huge resistivities and lack of semi- 
classical mean free path t (the putative mean free path 
would be smaller than the lattice spacing I < a; never- 
theless such "intrinsic" quantum diffusion can still be de- 
scribed by a diffusion constant extracted from the Kubo 
formula*™), whereby semiclassical description and per- 
turbative methods, based on the expansion in a small 
parameter l/kp£, break down; and (c) the Anderson 
localized regime when disorder becomes strong enough 
to push the conductance of a disordered sample below 4 
2e 2 /h. Note that to observe the effects stemming from 
localization of wave functions, the size of the conductor 
has to be greater than the localization length L 3> £ — 
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on length scales smaller than £ one cannot differentiate 
an Anderson insulator from a disordered metal. Fig- 
ure delineates the boundaries of these regimes for a 
system modeled by the half-filled Anderson Hamiltonian 
of Eq. (UJ. 

The transmission matrix t 



t = 2V -ImEz, • G\ 



N 



-Im E 



R, 



is obtained from the real-space Green function G r ' a 

1 



G r ' a = 



E -H - Y, r > a 



(2) 



(3) 



where G\ N , G a N1 are submatrices of G r - a {G a = [G r ] t ) 
that connect layers L = 1 and L = N of the sample 
along the x-axis. Here ImEi^ = (££ R — T,1 R )/2i are 
self-energy matrices (r-retarded, a-advanced) which de- 
scribe the coupling of the sample to the leads^ with E r = 
t r L + t r R (E a = [E r ]t). This particular computationally 
efficient implementation of the Landauer-Buttker formal- 
ism, which takes the microscopic Hamiltonian as an in- 
put, has its origins in the treatment of tunneling current 
in metal/insulator/metal (MIM) junctions — it was de- 
veloped in order to evade pathological properties of a 
tunneling Hamiltonian when attempting to take into ac- 
count higher order tunneling processes. 17 

All of the results shown in Sec. 11111 and Sec. II VI are ob- 
tained by evaluating exactly the Landauer transmission 
matrix for zero-temperature quantum transport in the 
half-filled (Ep — 0) Anderson Hamiltonian Eq. Q for a 
finite-size barrier. The disorder averaging is performed 
over an ensemble containing 1000 different samples for 
metallic disorder strengths W < 16.5 and, due to the 
need to search for rare events T n ~ 1 in special configu- 
rations of disorder, for 10000 samples on the insulating 
side W > 16.5. 



III. TRANSMISSION THROUGH DISORDERED 
INTERFACES AND THIN BARRIERS 

The distribution function of the eigenvalues T n is for- 
mally defined as 



P(T) = (J28(T-T n 



(4) 



where (...) stands for averaging over all possible real- 
izations of impurity configurations for a given disorder 
strength. Early mesoscopic studies of phase-coherent dis- 
ordered conductors have been focused on bulk systems in 
the weak scattering regime, where one finds celebrated 
perturbative quantum interference effects (such as weak 
localization and conductance fluctuations) within diffu- 
sive transport regimei For such systems, an analytical 
expression for P(T) has been obtained for the first time 
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FIG. 2: The distribution of transmission eigenvalues P(T) 
obtained in an ensemble of 1000 disordered barriers, at each 
junction thicknesses L. The chosen disorder strengths W of 
the random potential in the Anderson Hamiltonian generate 
the following systems in thick enough barriers (Fig. [^: (a) 
for W = 4, diffusive semiclassical metal (b) for W = 10, bad 
metal (c) for W = 22, Anderson insulator (here we use an 
ensemble of 10000 barriers). The dashed and dotted line plot 
the Schep-Bauer Psb(T) and the Dorokhov Pu(T) distribu- 
tions, expected to be valid for dirty interface and diffusive 
semiclassical metal, respectively. Note that these are not fits, 
but analytical expressions [see Eq. JSJ and Eq. JBJ] that de- 
pend on the disorder-average barrier conductance as a single 
parameter. 
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FIG. 3: The distribution of transmission eigenvalues P(T) 
and the conductance (at Ef = 0) of tunnel barriers of dif- 
ferent thickness that are made of disordered binary alloy 
A0.5B0.5 [panel (a)]. The binary alloy is modeled by the 
Anderson Hamiltonian Eq. Q with ea = ~£b = 4.6 being 
randomly distributed on its diagonal. This random poten- 
tial energy induces a hard gap in the eigenspectrum around 
Ef = 0, as shown in the panel (b), which is phenomenolog- 
ically similar to the gap in the density of states of band or 
Mott insulators. 



by Dorokhcn*i& 



Pd(X) 



(G) 



g q tVT^t" 



(5) 



and rederived within different theoretical frameworks 
Here (G) is the disorder-averaged conductance. The dis- 
tribution Prj (T) is universal in the sense that it does not 
depend on sample-specific properties (such as dimension, 
geometry, and carrier-density) . Although strictly derived 
for a quasi-one-dimensional wire (i.e., wire whose length 
is much bigger than its width), the scaling of transmis- 
sions implied by Pd (T) seems to have much wider valid- 
ity, as long as the conductor is in the (Bloch-Boltzmann) 
metallic regime^ 

The importance of interface scattering in giant magne- 
toresistance phenomenal has given an impetus to reex- 
amine transport through disordered interfaces. For the 
transparency of dirty interface, whose disorder-averaged 
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FIG. 4: The number of open conducting channels, whose 
transmission eigenvalues are close to one T n 6 [0.95, 1], in 
an ensemble of 10000 barriers for a given junction thickness. 
The barriers are made of strongly disordered materials, char- 
acterized by the disorder strength W = {16.5, 22, 30}, which 
behaves as an Anderson insulator in the bulk (see also related 
Fig. [T]and Fig. HJ. 



two-probe conductance is much smaller than the conduc- 
tance of corresponding point contact (G) < Gqpc, a 
Schep-Bauer distributio n 21 ! 22 has been found to be ap- 
plicable 



Psb(T) 



(G) 



1 



(6) 



While Psb(P) has been derived 2 ^ in the limit where 
barrier thickness is much smaller than the Fermi wave- 
length Xp, recent experiments^ have suggested that it 
might be valid even for thicker strongly disordered bar- 
rier Xp <C L < £, on the proviso that its width is smaller 
than the localization length £. 

Besides diffusive wires and dirty interfaces, analytical 
expressions for P(T) has been found for chaotic cavi- 
ties and double barrier junctions, as well as for combi- 
nations of these four generic cases^ Despite important 
insights obtained from different approaches^i^iiSiSiiS 
that yield Pd (T) and Psb (T) , no theory exist that would 
make it possible to obtain explicit expression for P(T) of 
a 3D mesoscopic disordered conductor which is in the 
non-semiclassical diffusive regime (i.e., the bad metal in 
Fig. extending all the way into the localized regime. 4 

We plot in Fig. [21 numerically exact P(T) as a func- 
tion of the barrier thickness, obtained by diagonalizing 
ttt in Eq. for each sample of an ensemble of disor- 
der configurations. For metallic diffusive barriers [panel 
(a)], £ -c L <C £ (note that a quasi-one-dimcnsional sys- 
tem will inevitably turn into an insulators when L > £ 
independently of the strength of the disorder^), trans- 
mission eigenvalue distributions follow Pd (T) prediction. 
However, for barriers made of the bad metal, P(T) is not 
equal to either Pd(T) or weak-localization-corrected24 
Pd(P), even though it remains bimodal distribution with 
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most of channels being either closed T n ~ or open 
T n ~ 1 [panel (b)]. Note that no single interface on 
the metallic side of the LD transition W < 16.5 can be 
described by Psb(T') [panels (a) and (b)]. 

When the disorder is strong enough to drive the LD 
transition in the bulk 3D samples, Psb{T) becomes valid 
in the single plane limit [panel (c)]. Moreover, it is also 
useful to some extent to describe P(T) for barriers com- 
posed of few such planes, as suggested by experiment and 
semi-intuitive arguments of Ref.Oj Finally, in Sec. II VI we 
demonstrate that the shot noise provides very sensitive 
tool to compare different distributions P(T) encountered 
here, as well as to differentiate those that are apparently 
similar [such as the distributions in panel (c)]. This is 
due to the fact that non-trivial features of P(T), such 
as the appearance of open channels in disordered tun- 
nel barrier^ii directly affect the suppression of the shot 
noise power below its trivial (Poisson limit) value which 
characterizes pure tunneling. 

To contrast the transport through strongly disordered 
barriers with tunneling through barriers made of a mate- 
rial with a gap in the density of states (such as the con- 
ventional band-insulators determined by single-particle 
quantum mechanics^ or more intricate Mott insulators 
which are governed by strongly correlated physics^ 5 .), we 
introduce disordered binary alloy A0.5B0.5 between the 
metallic leads. This system, which is composed of an 
equal number of atoms A and B randomly distributed 
throughout the simple cubic lattice, is modeled by ran- 
dom potential energy ea = — £b on the diagonal of TBH 
in Eq. (JTJ. As shown in Fig. |3 large enough ea — \sb\ 
will open a hard gap around Ep = in the density of 
states (DOS). The single interface of a solid with the gap 
in the bulk DOS can display distinctive transport prop- 
erties^ which manifest here as a non-trivial distribution 
P(T) where non-negligible transmission eigenvalues have 
finite probability to appear. However, already for the ul- 
trathin barriers L — 4, all transmission eigenvalues fall 
within the interval T n £ [0, 0.004], while the conductance 
exhibits typical exponential decay as a function of L. 

On the other hand, the disordered barriers always dis- 
play a non-trivial distribution of transmission eigenval- 
ues, which can accommodate open channels even at very 
large W and beyond the ultrathin limit. In the case of 
single interfaces and ultrathin barriers, the disordered 
region does not provide enough spatial extension in the 
transport direction to allow for the localization of wave 
functions i2i This also leads to emergence of the low- 
energy extended states within conventional ultrathin alu- 
minum oxide barriers that contain disorder or defects. 5 
We plot in Fig. 0] the decay of the number of open chan- 
nels as the barrier thickness increases, where the disorder 
strengths correspond to the Anderson insulator in Fig.^ 
The appearance of open channels beyond ultrathin bar- 
rier widths is a type of a rare event in the Anderson in- 
sulating phase (note that other types of rare events can 
arise in special configurations of disorder, even in the 
metallic phased). 



When Anderson insulator samples become larger than 
L^, phonon-assisted tunneling allows charges to propa- 
gate by hopping between the localized sites thereby gen- 
crating a finite conductance. However, the transport 
studied here takes place through phase-coherent barriers 
(i.e., their size satisfies L, L V ,L Z < L^) and, therefore, ef- 
fectively at zero temperature. The open channels inside 
the Anderson insulator junctions are due to the tunneling 
via rather special configurations of localized states that 
provide a path for resonant transmission of electrons. 29 
One example of such rare event is a wave function, with 
energy close to the Fermi energy (Ep = 0), which is 
symmetric with respect to the leads. Such wave function 
would make possible resonant transmission T n ~ 1, so 
that the conductance is proportional to the probability 
of finding such special barrier. This can be seen by com- 
paring Fig. 0] (which essentially gives the probability to 
encounter an open channel in a given ensemble of barri- 
ers) to the corresponding barrier conductances plotted in 
Fig.0 

IV. LINEAR STATISTICS: SHOT NOISE AND 
PROXIMITY CONDUCTANCE 

Over the past decade experimental and theoretical 
investigation of the shot noise, as a random process 
characterizing non-equilibrium state into which a phase- 
coherent conductor is driven by the applied voltage, has 
become one of the most active frontiers in mesoscopic 
physics^ The power spectrum of the shot noise, at zero 
frequency and at zero temperature, can be expressed^^ 
in terms of the Landauer transmission eigenvalues T n for 
non-interacting electrons transported through a conduc- 
tor attached to two leads 

r+°° _„ o P 2 £ " x f^ 

S = 2 dt' [I(t)I(t')-f] = 2eV^- T n (l-T n ). 

" / -°° h n=l 

(7) 

Here I is the time-average of the current flowing through 
the system under the applied voltage V . Thus, by mea- 
suring the shot noise one effectively probes second mo- 
ment of P{T), thereby obtaining complementary infor- 
mation to traditional conductance that is associated with 
the first moment of P(T). 

The suppression of the shot noise power S — 2FeI 
with respect to the Poisson limit S — 2el is quanti- 
fied by the Fano factor F. In the pure tunneling regime 
T n <C 1 => F = 1 because transfer of electrons through 
the barrier is uncorrelated in time and, therefore, de- 
scribed by the Poisson statistics. On the other hand, in 
the diffusive metallic conductors (more precisely, in the 
disordered Bloch-Boltzmann conductors in Fig. ^ whose 
size is such that £ <C L), the shot noise power is re- 
duced by a factor F = 1/3. This is due to the corre- 
lations generated by Fermi statistics — electron injection 
into the conductor is less likely if another electron is al- 
ready occupying one of the conducting channels. In the 
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FIG. 5: The disorder- averaged Fano factor, quantifying sup- 
pression of the shot noise power S = 2FeI from its Poisson 
value F = 1, as a function of the barrier thickness. Each 
curve is parametrized by the strength of the disorder W intro- 
duced in barrier (see Fig. . Note that the three horizontal 
lines label: (i) F = 1/3 shot noise suppression expected in 
the diffusive semiclassical conductors [i.e., Bloch-Boltzmann 
metal in Fig. Q whose transparency is described by Pu(T)]; 
(ii) F = 1/2 for dirty interfaces described by Psb(T); and (iii) 
F — 1 as a signature of pure tunneling through an insulator. 
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ballistic limit T n = 1, Pauli principle correlating non- 
interacting fermions leads to a complete noise suppres- 
sion F = 0. The shot noise is a genuine quantum trans- 
port phenomenon since deterministic classical transport 
also suppresses S to zero due to the lack of stochasticity 
associated with quantum mechanical propagation of elec- 
trons. The "magic" suppression factors, such as F — 1/3, 
are expected to be universally valid, i.e., independent on 
the details of the system such as geometric parameters 
of the conductor or its resistance. Since F = 1/3 follows 
from Pd(T) used in Eq. Q, while P S b(T) gives F = 1/2, 
the Fano factors may serve as an indirect and experimen- 
tally observable confirmation of a particular distribution 
P(T). 

The knowledge of P(T) makes it possible to compute 
the disorder-average of any quantity that can be cast into 
a form of the so-called linear statistics A — J2 n a (^n) 



(A) = r£a(T n )\ = J dT a(T)P(T). 



(8) 



The most frequently investigated examples of such quan- 
tities^ measured in the two-probe geometry and at zero- 
temperature, are: 

(a) the Landauer conductance 



<G> 



(b) the Fano factor 

(EnW-Tn)) _ fdTT(l-T)P(T) 



(F) = 



J dTTP(T) 



(9) 



(10) 



FIG. 6: The full distribution function of the Fano factor for 
several ensembles of barriers whose disorder-averaged (F) is 
plotted in Fig. 



(c) the linear conductance of 
region/superconductor (NS) junction 



a normal- 



(Gns) = 




P(T), (11) 



which holds in the zero-voltage, zero-temperature, and 
zero-magnetic-ficld limit, and for disorder confined to the 
N region. When the transparency of the NS interface is 
small (e.g., due to an insulator in between), single parti- 
cle tunneling is the dominant transport mechanism which 
renders (Gns)/(G) <C 1 ((G) is the conductance of the 
junction in the normal state). This is due to the fact 
that there are no available states within the energy gap 
A of S. However, in disordered- metal/superconductor 
junctions with transparent NS interface, (Gjvs) is en_ 
hanced due to the proximity effect which is microscop- 
ically generated by Andreev reflection at the NS inter- 
face. In this process, an incident electron is reflected as a 
hole, while a Cooper pair is pushed into the superconduc- 
tor. The expression for (Gns) is obtained by taking into 
account Andreev processes via Bogoliubov-De Genncs 
equations, while neglecting the self-consistency issues^ 
(e.g., superconducting order parameter is assumed to be 
a step function, thereby neglecting its depression on the 
S side of the junction^ as well as the terms of the order 
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FIG. 7: The disorder-averaged (over 10000 samples) conduc- 
tance of the Anderson insulator junctions of different thick- 
ness, attached to two metallic leads [panel (a)]. The panel 
(b) plots the ratio (Gns) /{G) of the disorder-averaged lin- 
ear conductance of a normal-region/superconductor junction 
[where normal-region is the barrier from panel (a)] and (G) 
from panel (a). The dotted horizontal line serves to high- 
light that (Gns) I (G) is mostly confined within the interval 
[0.4,0.6]. 



(A/ Ep) 2 ). For non-interacting quasiparticles that par- 
ticipate in purely Andreev processes at a perfectly trans- 
parent NS interface Gns/G < 2 [the upper bound is set 
by the ratio of Eq. ijTTJl and Eq. for T n — 1] . Another 
superconducting technique which allows to experimen- 
tally probe the transparency of atomic^ and mesoscopic 
conductors^ is to sandwich them between two supercon- 
ducting leads and analyze the subharmonic gap structure 
of the I — V characteristic of such Josephson junctions, 
which turns out to be determined" by P(T). 

Figure |3] demonstrates that F = 1/3 suppression 
is indeed applicable in the Bloch-Boltzmann transport 
regime, i.e., in the transport through thick enough bar- 
riers (but still L < L y ,L z ) where semiclassical diffusive 
(£ <C L) charge propagation takes place. Moreover, the 
suppression factor saturates (F(L)) — > Fm as a func- 
tion of the barrier thickness also for diffusion through the 
bad metal barrier. However, its asymptotic value Fm is 
steadily increasing Fm > 1/3 as a function of W when 



G NS /G 



FIG. 8: The full distribution function (sampled over 10000 
realizations of disorder) of Gns/G (see also Fig. [TJ quantify- 
ing mesoscopic sample-to-sample fluctuations for the Ander- 
son insulator barriers characterized by the disorder strength 
W = {22, 30} and thickness L = {1, 6, 14}. 



W > 6 crosses over the boundary of semiclassical trans- 
port regime in Fig. ^ For the Anderson insulator barri- 
ers W > 16.5, the Fano factor is 1/3 < (F(L)) < 1, as 
long as there is a probability to encounter open channels 
(Fig. through the barrier. It is also monotonic func- 
tion of the barrier thickness since P(T) scales with L in 
a fashion shown in panel (c) of Fig. [5] When all channels 
become closed in thick Anderson insulator barriers, the 
Fano factor reaches its trivial asymptotic value (F) = 1, 
thereby signaling that pure tunneling takes place through 
such barriers. Thus, Fig. suggests that the Fano factor 
offers a unique single scalar quantity that is able to re- 
solve disordered thin barriers with different transmission 
properties, as well as to label diffusive transport regimes 
of Fig. ^ within thick barriers. 

Quantum coherence, its non-local features, and ran- 
domness of microscopic details cause large fluctuations 
of physical quantities in disordered mesoscopic systems^ 
Contrary to the intuition developed from thermal fluctu- 
ations (and their self-averaging properties) in statistical 
physics of macroscopic systems, the average value and 
variance are not enough to characterize the distributions 
of various physical quantities in open (e.g., conductance, 
local density of states, current relaxation times, etc.) 
or closed (e.g., eigenfunction amplitudes, polarizability, 
level curvatures, etc.) mesoscopic systems. These distri- 
butions can become particularly broad upon approaching 
the LD transition^ Thus, Fig.EJintroduces the full distri- 
bution function of the Fano factor revealing that even in 
the diffusive metallic barriers there are sample-to-sample 
fluctuations yielding wide distributions (over the interval 
F e [0, 1]) when disorder is increased. 
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The possibility of fully open channels (with trans- 
parencies close to one T n ~ 1) to appear in the An- 
derson insulator barriers, as demonstrated in Sec. II I II 
has been indirectly suggested by a zero-bias anomaly 
in the I — V characteristic of normal- metal/ Anderson- 
insulator/superconductor (NTS) junctions^ Further- 
more, the comparison of the conductance (Gns) of the 
NIS junction with the conductance (G) of the junction 
in the normal state (i.e., (G) is the conductance of the 
metal/ Anderson-insulator/metal junction) allows one to 
test the importance of different special configurations of 
localized states that make possible resonant tunneling re- 
sponsible for the increase of the zero-bias conductance. 
For example, in the case of resonant tunneling through 
a chain of two localized sites^i (Gns) ~ 0.27(G). The 
same trend was conjectured to persist in barriers where 
tunneling through quasi-one-dimensional chains of arbi- 
trary number of localized states can occur i2i As the thick- 
ness of the barrier increases, more complex configura- 
tions would allow for percolation paths through localized 
states. However, they have not been observed in exper- 
iments measuring the normal conductance (G). Never- 
theless, the puzzling finding of experiments^ on NIS 
junctions is Gns/G G [1.0, 5.0]. 

Figure plots the conductance of the Anderson in- 
sulator junctions, as well as the corresponding ratio 
{Gns) /(G) when one of the normal leads of the MIM 
junction is turned into a superconducting one. The distri- 
butions of P(T) obtained in Sec. ITTT1 yields (G NS )/(G) G 
[0.4, 0.6] as a function of the barrier thickness. However, 
we recall here that conductance fluctuations in strongly 
disordered phase-coherent samples can reach the same 
magnitude as the conductance itself4*22 Therefore, we 
investigate full distribution function of Gns /G in Fig. [SI 
which shows that particular phase-coherent samples can 
indeed exhibit 0.27 < Gms/G < 2 similarly to the ones 
found in experiments^ Nonetheless, on many NIS junc- 
tions of Ref. ^ Gns/G > 2 is observed. This sug- 
gests that the interplay of proximity effect in the An- 
derson insulator and electronic interactions^ 3 , (that can 
play an important role in the localized phase due to lack 
of screening) takes place. Such effects are not captured 
by Eq. Hll|) that takes into account only Andreev reflec- 
tion of non-interacting quasiparticles at the NS interface. 
Their treatment would require more involved theoretical 
approaches, such as possible combination of dynamical 
mean-field theory extended to inhomogeneous systems^ 3 - 



(that include superconducting regions)^ 3 - with the typical 
medium theory of Anderson localization 3 ^ which would 
make it possible to study proximity effect in strongly cor- 
related and disordered systems (modeled by the standard 
Hubbard model with diagonal disorder used here). 



V. CONCLUSION 

We have investigated how statistics of the Landauer 
transmission eigenvalues P(T) for 3D barriers attached 
to two ideal metallic leads scales with the thickness of 
the barrier, as well as its dependence on the disorder 
strength which determines different quantum-transport 
regimes. When barriers are made of the bad metal (char- 
acterized by exceptionally high resistivity and lack of 
semiclassical mean free path), P(T) remains bimodal, 
but it does not obey scaling predicted by the standard 
Dorokhov distribution. The validity of the Dorokhov dis- 
tribution is confirmed for conductors where semiclassical 
diffusive metallic transport takes place, but which are not 
just quasi-one-dimensional wires of length much greater 
than its cross section, as assumed in different theoretical 
derivations. The characteristic signature of the distri- 
butions of all metallic (semiclassical or quantum) diffu- 
sive barriers is encoded into the scale independent Fano 
factor F > 1/3 measuring suppression of the shot noise 
power. In special configurations of disorder, strongly dis- 
ordered (i.e., Anderson insulator) barriers can accommo- 
date fully open channels T n ~ 1 due to resonant trajecto- 
ries through localized states. In experiments, this would 
lead to systems such as MIM junctions with F < 1 or 
normal- metal/ Anderson- insulator/superconductor junc- 
tions where ratio Gns/G takes any value < Gns/G < 
2 allowed within the proximity theory that excludes elec- 
tron correlation effects in the Anderson insulator phase. 
On the other hand, the explanation of Gns/G > 2 would 
require to treat proximity effect in strongly correlated 
and strongly disordered systems. 
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